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Abstract 

We study the tunneling through delta and double delta po- 
tentials in fractional quantum mechanics. After solving the 
fractional Schrodinger equation for these potentials, we cal- 
culate the corresponding reflection and transmission coeffi- 
cients. These coefficients have a very interesting behaviour. 
In particular, we can have zero energy tunneling when the 
order of the Riesz fractional derivative is different from 2. 
For both potentials, the zero energy limit of the transmis- 
sion coefficient is given by 7o = cos 2 (it /a), where a is the 
order of the derivative (1 < a < 2). 



1. Introduction 

In recent years the study of fractional integrodifferential equations applied to physics 
and other areas has grown. Some examples are [TJ |2j [3], among many others. More 
recently, the fractional generalized Langevin equation is proposed to discuss the anoma- 
lous diffusive behavior of a harmonic oscillator driven by a two-parameter Mittag-Leffler 
noise [I]. 

Fractional Quantum Mechanics (FQM) is the theory of quantum mechanics based 
on the fractional Schrodinger equation (FSE). In this paper we consider the FSE as 
introduced by Laskin in [5J |B]. It was obtained in the context of the path integral 
approach to quantum mechanics. In this approach, path integrals are defined over Levy 
flight paths, which is a natural generalization of the Brownian motion [7] . 

There are some papers in the literature studying solutions of FSE. Some examples 
are [H [91 [10]. However, recently Jcng ct al. [llj have shown that some claims to 
solve the FSE have not taken into account the fact that the fractional derivation is a 
nonlocal operation. As a consequence, all those attempts based on local approaches are 
intrinsically wrong. Jeng et al. pointed out that the only correct one they found is the 
one [12] involving the delta potential. However, in [12] the FSE with delta potential was 
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studied only in the case of negative energies. This has been generalized in [13] , where we 
have solved the FSE for the delta and double delta potentials for positive and negative 
energies. 

The objective of this paper is to study the tunneling through delta and double delta 
potentials in the context of the FSE. As a result, we found some very interesting prop- 
erties that are not observed in the usual a = 2 quantum mechanics. Probably the most 
interesting is the presence of tunneling through delta and double delta potentials even 
at zero energy. Moreover, in the case of the double delta potential, this zero energy 
tunneling is independent of the relation of the two delta functions. In Lin et al. [2] the 
problem of calculating the transmission coefficient in FQM for the double delta potential 
has been addressed; however, the authors have used that same local approach that Jeng 
et al. [TTJ have shown to be wrong. As expected in this case, our results differs from 
theirs. 

We organized this paper as follows. Firstly, and for the sake of completeness, we 
reproduce the solution of the FSE for the delta and double delta potentials, as given in 
|13j . presenting their respective solutions in terms of Fox's if -function. Some calculations 
and properties of the Fox's if-function are given in the Appendixes. Then we study 
the asymptotic behaviour of those solutions, calculate the reflection and transmission 
coefficients, and study some of their properties. The limit a — > 2 for these coefficients 
and the boundary conditions satisfied by the solutions are also discussed in two other 
Appendixes. 

2. The Fractional Schrodinger Equation 

The one-dimensional FSE is 



di/;(x, t) 



D a (-h 2 A) a/2 il)(x,t) + V(x)^(x,t), 



(1) 



where 1 < a < 2, D a is a constant, A = d 2 is the Laplacian, and (— h 2 A) a ^ 2 is the Riesz 
fractional derivative [12], that is, 





(3) 



The time-independent FSE is 



D a (-h 2 /\) a/2 ip(x) + V(x)i/>(x) = E^{x). 



(4) 



In the momentum representation, this equation is written as 



(5) 
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where (W * <f))(p) is the convolution 

/+oo 
W(p-q)<P(q)dq, (6) 
-oo 

and W(p) = J-"[V(x)] is the Fourier transform of the potential V{x). 

Solutions of the FSE for delta and double delta potentials are given in [T3] in the 
situations of bound and scattering states. Since we need to study the asymptotic be- 
haviour of these solutions in order to find the transmission coefficients, and for the sake 
of completeness, we will reproduce the calculations of the wave functions in the case of 
scattering states. 



2.1. FSE for Delta Potential 

Let us consider the case 

V{x) = V S(x), (7) 

where 5(x) is the Dirac delta function and Vq is a constant. Its Fourier transform is 
W(p) = Vq and the convolution (W * <p){p) is 

(W*<j>)(p) = V K, (8) 

where the constant K is 

/+oo 
<P(q)dq. (9) 
-oo 

The FSE in the momentum representation §5§ is 

\p\ a - Hp) = -iK, (io) 

where 

Since we are interested in scattering states, we will consider that E > and write 

7^ = A Q , (12) 
where A > 0. Since f(x)S(x) = f(0)5(x), the solution of Eq. fnoj) in this case is 

*W = Ip^-X* + 2lThC ^(P - A ) + ^C 2 5(p + A), (13) 

where C\ and C 2 are arbitrary constants and the constant 27rh was introduced for later 
convenience. Using this in Eq.(|9]) gives that 

K = ->yK J +C ° Xa + toWx + 2nhC 2) (14) 
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where the integral is interpreted in the sense of Cauchy principal value, and it gives 
,+ °° dp 



\p\ a - A c 



2A 



l-a 







dq 
q a - 1 



7T 7T 



-2A^ Q -cot-, 

a a 



(15) 



where we have used formula 3.241.3 (pg. 322) of [16J - see Eq. fl89|) . The constant K is 
therefore 

2nh(C 1 + C 2 )a\ a - 1 

(16) 



K 



and we have 

<j){p) = 27ihC 1 5(p - A) + 2nhC 2 6(p + A) 



a\ a 1 — 27T7 cot (7r/a) 

27rfr 7 (Ci + C 2 )aA Q - 1 



(17) 



(aA"- 1 - 2tt7 cot (vr/a)) (|p| Q - X a ) 
Next we need to calculate the inverse Fourier transform of <p(p) to obtain ip(x), that is 



^( x ) = de iXx/n + C 2 e~ iXx/h 



2it 1 a(C l + C 2 ) fXx\ 
(aA«- 1 -27T7C0t(7r/a)) a \~h ) 



where 3 a ( w ) is the Cauchy principal value of the integral 

i p+oo 



and such that 



+oo ipx/h 



cos wq 



n Jo <t 



dq, 



dp = 2ir\ L - a 3 a (\x/h). 



(19) 



(20) 



, \p\ a — A c 

The above integral is calculated in the Appendix B, and is given by Eq.( 19~Tj) . Using this 
result, and the definition of 7 in Eq. fllip and A in Eq. (jl2j) . we can write that 



(21) 



where 



\\x\ 



ah , H 2 t i 



X\x\ 



'2,3 



h 



-H, 



2,1 
2,3 



X\x\ 



\\x\Y (l,l),(l,(2 + a)/2) 

(l,a),(l,l),(l,(2 + «)/2) 

(l,l),(l,(2-a)/2) 
(l,a),(l,l),(l,(2-«)/2) 



(22) 



with H. 



2.1 
2.3 



denoting a Fox's if-function (see Appendix A), and 



a 



E\ ~ 7T 

— - cot - 

U a 



and 



ahD 



l/a 



a/(a-l) 



(23) 



(24) 



4 



2.2. FSE for Double Delta Potential 

Now let the potential be given by 

V(x) = V [5(x + R/2) + /j5(x - R/2)], (25) 

with /i, Vq and R real constants. When Vq < this potential can be seen as a model for 
the one-dimensional limit of the molecular ion H 2 [T7] . The parameter R is interpreted 
as the internuclear distance and the coupling parameters are Vb and /iVo. Its Fourier 
transform is 

W(p) = V e ipR/2h + V ne- ipR/2R (26) 

and for the convolution 

(W*<f>)(p) = V^ pRj2h K x (R) + V 0f ie- ipR/2h K 2 (R) , (27) 
where Ki(R) and K 2 (R) are constants given by 

/+oo 
e- iRg/2h <f){q)dq. (28) 
■oo 

The FSE in momentum space is 

" J~) = ~^ Rp/2hK ^ R ) ~ l^- tRp/2R K 2 (R), (29) 

where we used the notation introduced in Eq. lfiTT) . 

Since we are interested in scattering states, we have E > and we write A as in 
Eq. (1121 and for the solution of Eq. (T291) we have 



m = ~ A) + 2«W 2 6( P + A) - 2!^W _ /fTg^W , 

v 7 v 7 v 7 - X a \p\ a - X a v 7 

Using this expression for <p(p) in Eq. (|28|) of definition of K\(R) and K 2 [R) we have 

(1 + 27r 7 A 1 ^L(0))iri(i?) + /i2vr 7 A 1 ^ a (Ai?/n)ir 2 (i?) = 2nhC[, (31) 
2vr7A 1 - a a a (A J R/n)K 1 ( J R) + (1 + ^t^A^^O))/^) = 2nhC' 2 , (32) 

where 

= C ie - lRX / 2h + C 2 e iRX / 2h , C 2 = C x e iRX ' 2h + C 2 e~ iRX l 2h . (33) 
In order to write the solution of these equations it is convenient to define 



27T7 a \U 

where U was defined in Eq. (I24[) . in such a way that we have 

K X {R) = ^£[(^-1+^(0))^ -Z a {\R/K)C' 2 \, (35) 
K 2 (R) = ^[( e + 3a(0))C£ -3 a (\R/h)C[}, (36) 
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where 

W = (e + X(0))(^" 1 + 3«(0)) - (X(Ai?/^)) 2 . 
Using Ki(R) and K 2 (R) in Eq. (l30l) gives Then, for V'(x), we have 

^( x ) = Cie iA:c/rt + C 2 e~ iXx/h 

+ 2^[(e + 3«(0))(T 2 -3a(Ai2/fi)(r i ]$ a ( A|X ~/ /21 ) , 
where we have expressed the result in terms of the function $ a defined in Eq.( l22|) . 



(37) 



(38) 



3. Calculation of the Transmission Coefficients 

In order to calculate the transmission coefficients, we need to know the asymptotic 
bahaviour of the solutions. The asymptotic behaviour of Fox's ff-function is given, if 
A > 0, by Eq.flH2D or Eq.flBU) according to A* > or A* = 0, respectively - see Eq.flEU). 
In Q a (\\x\/h) we have the difference between two Fox's if-functions of the form 



H. 



2,1 
2,3 



W 



;i,a),(l,l),(l, yU ) 



for /i = (2 + a)/2 and — (2 — a)/2. In both cases we have A = a > 0, but A* = for 
/i = (2 + a)/2 and A* > for fi = (2-a)/2. Therefore, using Eq. QB2J when p = (2-a)/2 
and Eq.( l84")) when /i = (2 + a)/2 we have, respectively, that 



H. 



2.1 
2.3 



W 



H. 



2.1 
2.3 



(l,l),(l,(2 + «)/2) 
(l,«),(l,l),(l,(2 + «)/2) 

(l,l),(l,(2-«)/2) 
(1, a), (1,1), (1, (2 -a)/2) 



2w • 

— sinro + o(l 



w 



\w\ — > oo, 



\w\ — > oo, 



and then 



A|x| 
~7T 



2 sin ■ 



Alxl 



+ o(\x\~ 1 ), 



\x\ — )■ oo. 



(39) 
(40) 

(41) 



3.1. Transmission Coefficient for the Delta Potential 

The behaviour of the wave function ip(x) given by Eq.f l2ip for x — > ±oo is therefore 
Mx) = C 1 e iXxlh + C 2 e- iXxlh ±nJC l + C 2 )sixi^ + o(x- l ) ) x ±oo. (42) 



or 



^(x) = Ae iAa;//t + Be- iXx/h + o (aT 1 ) , x -oo, 
^(s) = Ce iAx//l + De"^ + o (x- 1 ) , x -> +oo, 



(43) 
(44) 
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where we defined 



A = C 1 + i(C 1 + C 2 )Sl a /2, B = C 2 - i{d + C 2 )Q a /2, (45) 
C = G x - z(Ci + C 2 )0 Q /2, D = C 2 + z(Ci + C 2 )a*/2. (46) 

Now let us consider the situation of particles coming from the left and scattered by the 
delta potential. In this case D = (no particles coming from the right) and B = rA and 
C = tA, where the reflexion 1Z and transmission T coefficients are given by 1Z = \r\ 2 
and T = \t\ 2 (see, for example, [IS])- The result is 



1 + z'Oq, 1 + i£l c 

and then 



(47) 



In Fig.(l) we show the behaviour of these coefficients for different values of a. This plot 
and the other ones of this paper have been done by means of numerical integration of 
Eq. fflUj) using Mathematica 7. 




Figure 1: Reflection and transmission coefficients as function of E/U, as given by 
Eq. pg]) . for different values of a. 

We must note that transmission coefficient has a very interesting behaviour at zero 
energy. If we take the limit E — > in the expression for Q a in Eq. f[2Uj) we see that 

7T 

lim Q a = — tan — , (49) 
e^o a 
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and then for the transmission coefficient T we have 



lim T = cos 2 — . (50) 
e^o a 

This is an unexpected and very interesting effect, which demands further interpretation 
(see Conclusions). 

3.2. Transmission Coefficient for the Double Delta Potential 

Let us introduce the following notations: 

u _ e^+ZajO) v _ e + 3a(0) x _ Za(XR/h) 

aW ' aW ' ' aW 1 ' 

The asymptotic behaviour of the wave function ip{ x ) given by Eq. fl38|) for x — > ±oo is 
therefore 



Mx) = C ie iXx / h + C 2 e- iXx ' h + (UCi - XC') sin X \ x + R I 2 \ 

n 

+ (VC 2 - XC[)sm X \ X ~ h R/2 ^ + o(aT 1 ), x^±oo, 



(52) 



or 



where we defined 



and 



ip{x) = A'e iXx/h + B'e~ iXx/h + o (aT 1 ) , x -»• -oo, (53) 
^[x) = C'e lXx/h + D'e~ lXx / h + o (aT 1 ) , x ->• +oo, (54) 



A' = Ci + M 1; B' = C 2 + M 2 , (55) 
C" = d - Afi, D' — C 2 — M 2 , (56) 



M x = z(pC 1 + aC 2 + zrC 2 ), (57) 
M 2 = -i{pC x + pC 2 - irC x ), (58) 



with 



W + V\ v XR (U + V\ XR v 
— X cos — — , a = — - — cos — X , 



2 J h V 2 / h 

U-V\ XR 
sin — — . 

n 



(59) 



As in the case of the delta potential, let us consider the situation of particles coming 
from the left and scattered by the double delta potential. In complete analogy we have 
D' — (no particles coming from the right) and B' = rA' and C = tA', where the 
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reflexion TZ and transmission T coefficients are given by K = |rf and T = \t\ 2 . The 
result is 



2(r + ia) (p 2 - a 2 - r 2 + 1) 



(p 2 — cr 2 — r 2 — 1) — 2ip' (p 2 — cr 2 — r 2 — 1) — 2ip 

and 7?. and T can be written as 



(60) 



v 4(a 2 + r 2 ) (p 2 - a 2 - r 2 + I) 2 

(p 2 - a 2 - r 2 + l) 2 + 4(a 2 + r 2 ) ' (p 2 - a 2 - r 2 + l) 2 + 4(<t 2 + r 2 ) ' 1 ' 

We can simplify these expressions a little bit once we note that 

_ T2= sinW) ; 
crVV 

t ^ita 2 ' (63) 

A 2 = l^_V±l3 (64 ) 

a (aW + sin^AiZ/fr)) 2 ' 1 ' 

In Fig. (2) we show the behaviour of the transmission coefficient for different values of 
a and in Fig. (3) for different values of p. 



and then 



with 




Figure 2: Transmission coefficients as function of XR/2H, as given by Eq. (l63l . for differ- 
ent values of a, when 27r7 =10 and p = 1. 
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Figure 3: Transmission coefficients as function of XR/2H, as given by Eq. (1631) . for differ- 
ent values of /i, when 27r7 = 20 and a = 1.8. 

As in the case of the delta potential, we have a very interesting behaviour for these 
coefficients when E 0. Firstly, using Eq. ([9"9~j) in the Appendix B, we have 

Z a (\R/h) = 3 a (0) + A l X a ~ 1 + A 2 X 2 + 0(A 3a_1 ), (65) 

where 

= (R/hr- 1 = {R/hf cot (37T / a) 

1 2r(a)cos(7ra/2)' 2 2a ' 1 ' 

Using this, the expression Eq. (l37|) for W gives 

W = ^(0)(^(l+^ 1 )-2A 1 )A Q - 1 + (i/V 1 -^ 2 )A 2( ^ 1) -2a Q (0)A 2 A 2 + O(A Q+1 ), (67) 
where H = l/2n r ). Then, with some calculations, we obtain that 

4a 2 W 2 (a 2 + r 2 ) = Bi\ 2{a ~ l) + B 2 X a+1 + 0(X 2a ) } (68) 

with 

B 1 = (H(1+ I 2~ 1 )-2A 1 ) 2 } B 2 = -2(H(1 + fi- 1 ) -2A 1 )(2A 2 +3 a (0)(R/fr) 2 ), (69) 
and 

a 2 W 2 {p 2 - a 2 - t 2 + 1) = ^A 2(Q - 1} + ^A 4(Q - 1} + B' 3 X a+l + 0(X 2a ), (70) 
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with 

B[ = a 2 3l(0)(H(l + fi- 1 ) - 2A 1 )\ B' 2 = a 2 {H 2 ^ 1 - A±) 2 , 
B' 3 = -2a3 a (0)(if(l + fx' 1 ) - 2A 1 )(2a 2 3 a {0)A 2 - {R/h) 2 ). 

Using these results we can easily see that, for E/D a = \ a — > 0, 



(71) 



lim A 2 a = 1 ( 72 ) 
e-+o a a 2 4 2 a (0) 

and, since 3a (0) = — (1/a) cot (ft/cx), that 

lim T = cos 2 — . (73) 

This the same result we obtained for the zero energy limit of the transmission coefficient 
for a single delta potential. Moreover, this limit does not depends on //, which is the 
parameter that relates the two delta functions in the potential in Eq. (l25jl . Again, this 
is a very interesting and unexpected result. 



4. Conclusions 



The tunneling effect in fractional quantum mechanics has some very interesting proper- 
ties which are not observed in the usual a = 2 quantum mechanics. The most interesting 
is the presence of tunneling through delta and double delta potentials even at zero en- 
ergy. Moreover, in the case of the double delta potential, this zero energy tunneling is 
independent of the relation of the two delta functions. Let us give a possible explanation 
to these results. 

FQM was defined from the point of view of path integrals. As well-known, in this 
approach, when the sum is taken over paths of Brownian motion type, we have the 
standard quantum mechanics. On the other hand, in FQM the sum is taken over the 
paths of Levy flights [5] , which are generalizations of Brownian motion, and such that the 
corresponding probability distribution has infinite variance. By means of Levy flights, 
there is a not negligible probability of a particle reaching faraway points in a single jump, 
in contrast to a random walk of Brownian motion type [19] . 

In FQM, an uncertainty principle still holds, but with an appropriated modification, 
that is, we have [5] 

(lAxH^dApr) 1 ^ > li<a> l <a < 2 , (74) 

and such that in the standard quantum mechanics limit we can also have /i = a = 
2. Thus, even when E — 0, the particle can have energy AE = (\Ap\v) 2 ^/2m and 
momentum (| Apl^) 1 ^. This may not be enough for tunneling through a delta potential 
in the standard case, but in the fractional one, where long jumps of Levy flights enter 
the sum in the path integral, this may be responsible for the tunneling with probability 
cos 2 (ir/a). 

Acknowledgements: We are grateful to the referees for many valuable suggestions. 
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A. Fox's //—Function 



The Fox's H— function, also known as if— function or Fox's function, was introduced in 
the literature as an integral of Mellin-Barnes type |20j . 

Let m, n, p and q be integer numbers. Consider the function 



Ylribi + Bis^rii-ai-Ais) 



Ms) 



i=i 



i=i 



Y[ F(l -bi-Bis) T( ai + A iS ) 



(75) 



i=m+l 



i=n+l 



with 1 < m < q and < n < p. The coefficients Ai and Bi are positive real numbers; cii 
and hi are complex parameters. 

The Fox's if— function, denoted by, 



H, 



p.q 



X 



y a p , A p ) 



(b q ,B q ) 

is defined as the inverse Mellin transform, i.e. 



H™> n (x) 



TTm,n 



X 



(Ol,AlJ 

(h,B^ 



, (ctp, Ap) 
,{KB q ) 



— lA(s)x-*ds 



(76) 



(77) 



where A(s) is given by Eq. fP75|) . and the contour L runs from L — zoo to L+ioo separating 
the poles of T(l — % — Ai s), {i = 1, . . . , n) from those of T(bi + Bis), {i = 1, . . . , m). 
The complex parameters a« and 6j are taken with the imposition that no poles in the 
integrand coincide. 

There are some interesting properties associated with the Fox's if— function. We 
consider here the following ones: 



P.l. Change the independent variable 

Let c be a positive constant. We have 



Tjm,n 

p,q 



x 



(a p , A p ) 
(b q ,B q ) 



cH m ' n 



x 



(a p , c A p ) 
ip q} cB q ) 



(78) 



To show this expression one introduce a change of variable s — )■ cs in the integral of 
inverse Mellin transform. 



P.2. Change the first argument 

Set aeM. Then we can write 



a Trm,n 

p,q 



X 



iflpi Ap) 

(b q ,B q ) 



TTm,n 

p,q 



x 



{cip -\- aAp, Ap) 
(b q + aB q , B q ) 



(79) 



To show this expression first we introduce the change a p — > a p + aA p and take s — > s — a 
in the integral of inverse Mellin transform. 
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P.3. Lowering of Order 

If the first factor (ax, Ai) is equal to the last one, (b q , B q ), we have 

(a 2 , A 2 ), • • • , {a p , A p ) 
(b l ,B 1 ),--- ,(6 H ,B H ) 

(80) 

To show this identity is sufficient to simplify the common arguments in the Mellin-Barnes 
integral. 



Tjm,n 



x 



(ax,Ax),-- 
(h, B 1 ) r -.,(b, 



dpi Ap) 

g-l, B q _ x 



){a 1 ,A 1 



TTm,n—l 
n p-l,q-l 



X 



P.4. Asymptotic Expansions 

The asymptotic expansions for Fox's iJ-functions have been studied in [21]. Let A and 
A* be defined as 

q p n p m q 

i=l i=l i=l i=n+l i=l i=m+l 

If A > and A* > we have 



H m n ( x ) = Y,[ h - x(ar ~ 1)/Ar + o(x^- 1)/Ar )], \x\^oo. (82) 



r=l 

where 



i Il?=i r(&j + (i - g r )B 3 /A r ) UU,^r r(i - a,- - (i - a,) VA-) 

and if A > and A* = we have 1221 



H^ n {x) = [h r x {ar ~ l)/Ar + o(x {ar ~ 1)/Ar )] 



r=l 



where 



+ Az> +1/2)/A (c exp[*(5 + Cx 1/A )] - d exp[-i{B + Cx 1/A )}) 

+ ( X (-+V2)/|A|) ; | x |^oO, 



c = (27ri) m+n - p exp [m( a r ~ X) 6 *)]' 

r=n+l i=l 
p m 

d = (-2m) m+n - p exp [ - 7T? ( a r - ^ fy) ttz] , 

r=n+l j=l 

i P 9 , / A A\ (H-l/2)/A 
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5 = n \a,\- a > n i^i b ', * = x> - 1% + 

1=1 j=\ j=\ j=l 



P.5. Series Expansion 



In J2U] we can see that in some cases there is a series expansion for Fox's if-function. 
For example, when the poles of YlJLi F(bj + Bjs) are simple, we can write 

m oo 

W-EE^ fe+,/)/Bi . (85) 

j=l u=0 

where 

(-1)- ngw, r(& t - BtiPj + u)/Bj) niU r(i - Oi + + ij/gj) 
jV y\Bi Ul m+ i r(i - + 5,(6, + ^/s,) nr =n+1 r(a, - ^ + • 1 j 



B. Calculation of the Integral in Eq.(H9|) 

Let Z a (w) be given by 

Taking the Mellin transform we have that 

i ttz r +o ° V~ Z 

M w {Za{w)]{z) = -T{z) cos — / ~^—dy. (88) 

n * Jo y ~ ± 

This last integral is given by formula 3.241.3 (pg.322) of [16J, that is, 

f + °° x^ 1 it an 

/ dx = -cot— , 89 

where the integration is understood as the Cauchy principal valu^l Therefore we have 

. , r~ / m / n 1 . . 7r(l — z) 7r(l — z) , , 

M w [Z a {w)]{z) = r z sin 1 - ; cot -i 1. (90) 

a 1 a 

Using the relation 2 sin A cos B = sin (A + B) +sin (A — B) and writing the sine function 
in terms of the product of gamma functions we can write that 

i r(*)r(i=*)r(i - 



M w [3 a (w)](; 



2ar((l-z)^i)r(l-(l-z)^) 

2ar((i-z)^)r(i-(i-z)^) 21 h 



3 We remember that in the inversion of the Fourier transform the integration is to be done in the sense 
of the Cauchy principal value [33] . 
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Taking the inverse Mellin transform and using the definition of the Fox's if -function we 
have that 



1 



2.1 



+ —#2,3 



2a 



w 



w 



(1 - 1/a, 1/a), (1 - (2 + a)/2a, (2 + a)/2a) 

(0, 1), (1 - 1/a, 1/a), (1 - (2 + a)/2a, (2 + a)/2a) 

(1 - 1/a, 1/a), (1 - (2 - a)/2a, (2 - a)/2a) 

(0, 1), (1 - 1/a, 1/a), (1 - (2 - a)/2a, (2 - a)/2a) 



+ 



Using the properties given by Eqs.( f78|79|) and replacing w by \w\ since $ a (— w) = Z a (w) 
we obtain 



2.1 



+ 



2\w 



■H. 



2,1 

2,3 



\W\ 



\W 



(l,l),(l,(2 + a)/2) 

(1, a), (1,1), (1, (2 + a)/a) 

(l,l),(l,(2-a)/2) 
(l,a),(l,l),(l,(2-a)/2) 



+ 



(91) 



Let us see what happens in the particular case a = 2. From the definition of Fox's 
if -function we can see that 



H. 



2.1 
2.3 



\W\ 



(1.1) , (1,0) 

(1.2) , (1,1), (1,0) 



(92) 



and that 



H. 



2.1 
2.3 



W 



(1.1) , (1,2) 

(1.2) , (1,1), (1,2) 



#1,2 



W 



(1,1) 

(1,1), (1,2) 



2 rj-1,1 



w 



(0,1) 

(0,1), (-1,2) 



But [20 



#1,2 



(0,1) 

(0,1), (1-6, a) 



Enb(z) 



(93) 



where E a ^(z) is the two-parameter Mittag-Leffler function. However, it is known [21] 
that 

E 2a(z) = 5i!^£. ( 94 ) 



Consequently, we have 



-"2,3 



W 



(1.1) , (1,2) 

(1.2) , (1,1), (1,2) 



Then for a = 2 we have 



and 



&2 



Xx 



M 2 #2,2(-M 2 ) 

\\x\ 



\w\ sm \w\. 



-sin , 
2 h 



+oo ^ipx/h 



\p\ 2 - A 2 



dp 



7T X\x\ 



(95) 

(96) 
(97) 
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We are also interested in the expression of 3 a {w) for small w. From Eq. fj9T|) we see 
that we need to know the behaviour of 



H. 



2.1 

2,3 



\W\ 



(1,1), (1,^) 
(l,«),(l,l),(l,/i) 



for small w. This is given by the series expansion from Eq. (l85|) . which gives 



H. 



2.1 
2.3 



1,1), (l,fi) 

r(i - l/q)r(l/a) z x i°- 

r(l — fi/a)T(fi/a) a T(l — 2/i/a)T{2fj J /a) a 
r(l - 3/a)r(3/a) , ^ 4/ 



r(l - 2/a)r(2/a) 2 



2/a 



(9* 



+ 



T(l — 3 / u/a)r(3 / u/a) 2a 



r(i-a)r(i) _ r(i-2a)r(2) 2 r(i-3a)r(3) z 3 4 
r(i - fj,)r(fj.) r(i - 2/i)r(2/i) r(i - 3^)r(3^) 2 1 

Using this in Eq.( l9ip we arrive, after some manipulations, to 

1 



3«H = 3a(0) + 



2T(a) cos (vra/2) 



„ , cot(37r/a;) 9 „. ^ 



2a 



where 



1 7T 
-Ja(O) = COt -. 

a a 



(99) 



(100) 



C. The Limit a = 2 



Let us calculate the transmission coefficients in the standard quantum mechanical limit 
and see that we recover the usual results. 

Firstly, let us consider the delta potential. The transmission coefficient is given by 
Eq. fl48p . so that we need to calculate i n this case, where Q a is given by Eq. (l23p . The 
result is that 

E 



a, 



-1 



u 



tuVq 
~2h^' 



(101) 



where we used the definition of U in Eq. f[24l) and D 2 = l/2m, and such that 

r = 



(mV£/2h 2 E) 



(102) 



which is the well-known result [IB]. 

Now let us consider the transmission coefficient for the double delta potential. Let us 
also consider the case /i = 1 since the result in this case is well-known [22]. In order to 
calculate A 2 in Eq. ([63l we need Eq. fl96l) . which by the way gives 



3 2 (0) = o. 



(103) 
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When fx — 1 we have 

« = V=^, t-0, (104) 

and then 

2 _ 6AW 2 a 2 

2 ~ [AW + sin 2 (XR/hW (1 ° 5) 
But in this case Eq.( l3"7|) together with Eq.f l96|) gives 

AW + sm 2 (XR/h) = Ae 2 . (106) 

From Eq.f l5ip and Eq. (l59|) we also have that 

9 / XR 1 XR\ 2 . _, 

AW 2 a 2 = e cos — + - sin — . (107) 

\ lb _ It/ 



Then for A 2 we have 



, .-, I 1 XR 1 XR \ 2 

A 2 = -cos — + — sin— . 108 



e 



ft 2e 2 ft 



Let us change the notation a little bit in order to compare with the standard result in 
the literature. Let us define 

» = w " = t (109) 

Using this notation in the above expression for A 2 , we have for T in Eq.( |63|) that 

r= — (no) 

(4 + \J3hq cos 2/i + (/3 2 /2) sin 2/x ] 2 ' v ; 

which is the result in [22], pg 160 (where /i = /i). 



D. The Boundary Conditions 

In [13] we have discussed how the Riesz fractional derivative can be written in terms of 
the Riesz potentials, that is, for < a < 1 we have 

(_A) a / 2 ^(x) = ^-1l l - a i){x), (111) 
ax 

and for 1 < a < 2 we have 

(-A)^(x) = -^-TZ 2 ~ a ij(x) } (112) 
dar 

where lZ a ip(x) is the Riesz potential of ^(x) of order a' given by [TB] 

= w x 1 , — rr / + °° , ^51 , de, (113) 

n ; 2r(«')cos(aV/2)y-oo k-ei 1 "" 1 ^ 
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for < a' < 1, and TZ a ip(x) its conjugated Riesz potential given by 

If we use these expressions for the Riesz fractional derivative in the FSE for the delta 
potential V(x) = Vq5(x) and integrate as usual from — e to +e and take the limit e — > 
we obtain that 

-ii>-y.(x) v ° 



o+ dx 



o- h a D t 



-m- (us) 



This condition and the continuity one ^(O - ) = ^>(0 + ) are the boundary conditions to be 
satisfied by i/j(x). 

There is an important point to be noted here: the expression 

^n 2 - a tfj(x) = ^-n^-Vipix) 

dx dx 

is not the Riesz fractional derivative of order a — 1, which is given, for < a — 1 < 1, 
by 

dx 

Therefore, it is wrong to write the condition ( |115p as 

(-A)^/^(x)\ 0+ - (-A) (a_1)/2 ^(x)| _ = ^V(O), (H6) 

with (— A) Q '/ 2 being the Riesz fractional derivative. Maybe in another version of FQM 
involving a fractional derivative defined in another sense, it can holds one such condition, 
but this is not the case when it comes to the Riesz fractional derivative. In [H] the 
authors have used a boundary condition of the above type in their attempt to solve the 
problem for the double delta potential. Besides the already discussed problem with their 
local approach, it also seems that the use of that boundary condition is not justified. 

Let us consider the solution in the case of the delta potential and show that it satisfy 
Eq.f lll5p . In [T3] we have seen that 

d i f°° 

—TZ 2 - a ij(x) = —— e- ipx/h \p\ a - 1 sign(p)(j)(p)dp. (117) 



dx y ' 2nh a j x 

Using 4>{p) given by Eq. (fT7]) we obtain, after calculating the integrals with the delta 
functions, that 

^ 2 ">(x) = 2id (j X 1 e lXx/h -2iC 2 (^Y V^/" 
dX W _ x W (118) 
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where 



(W 



„e*-l 



9 



-dg. 



(119) 



This integral can be calculated in a way analogous as we did in Appendix B. The result 



is 



7T 



sienfuO ( if. 



r2,l 
'2,3 



-H; 



2,1 
2.3 



(0,l),(0,(2 + a)/2) 
(0, a), (0,1), (0,(2 + a) /2) 

(0,l),(0,(2-a)/2) 
(0, a), (0,1), (0,(2- a) /2) 



(120) 



In order to study the limit x — > O 1 * 1 we need to use 

lim Hli 



\W\ 



(0,l),(0,/3) 
(0, a), (0,1), (0,0) 



a: 



(121) 



which can be calculated using the definition of the Fox's if -function and the residue 
theorem (the RHS comes from the residue at w = 0), and which gives 



Using this result, we have 



lim E a (w) =±£. 



a 



(122) 



Inn — 7?. 2 -'r(.r) ■ ±-(C\ ■ -C' 2 )fi„ [ j; 



a-l 



x^o ± dx ' 2 

On the other hand, ip(0) can be calculated using Eq.f llOOj) in Eq.f JTB"]) . which gives 

a-l 



(123) 



(d + C 2 )otQ a X a - 1 n a D a 
= ^ = t ~p oc(Cx + C 2 )il a 



(124) 



2vr 7 Vb 

where we used the definition of 7 in Eq.f fTT]) . We see, therefore, that Eq. (1115p is satisfied. 
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